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Abstract 

Stylized facts of empirical assets log-returns Z include the existence of (semi) heavy tailed 
distributions fz{z) and a non- linear spectrum of Hurst exponents r(/3). Empirical data con- 
sidered are daily prices of 10 large indices from 01/01/1990 to 12/31/2004. We propose a 
stylized model of price dynamics which is driven by expectations. The model is a multi- 
plicative random process with a stochastic, state-dependent growth rate which establishes a 
negative feedback component in the price dynamics. This Q-order model implies that the 
distribution of log-returns is Laplacian fz{z) ~ exp(— whose single parameter a can be 
regarded as a measure for the long-time averaged liquidity in the respective market. A com- 
parison with the (more general) WeibuU distribution shows that empirical daily log returns 
are close to being Laplacian distributed. The spectra of Hurst exponents of both, empirical 
data Ternp SLiid simulated data due to our model Ttheor, are compared. Due to the finding of 
non-linear Hurst spectra, the Renyi entropy (RE) R0{fz)ia considered. An explicit functional 
form of the RE for an exponential distribution is derived. Theoretical REs of simulated asset 
return trails are in good agreement with the RE estimated from empirical returns. 



*Contact address: sreimann@iew.unizh.ch 
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We regard a financial market as a large complex system. Then Statistical Physics implies that 
(some of its) macro observables, such as prices, might be independent from its micro realizations 
and therefore are common to almost all realizations of the system. Indeed the statistical behavior 
of price fluctuations on different markets exhibit much structure that is common to 'all' finan- 
cial markets. These properties are called 'stylized facts', for an enlightening survey see 0] and 
also the monographs T, "H". In turn, the existence of stylized facts suggests that price trails of dif- 
ferent assets are realizations of a more general random/complex system, called 'a financial market'. 

Prices are macro observables of a financial market in that their time evolution is generated 
by the successive trading activities of a huge amount of financial agents. This justifies to take a 
macroscopic perspective for modeling. The idea of our model is the following: 

People go the financial market to make their money work. They do so by investing their money 
into assets. If the agent has capital m to invest in asset A, he can buy |^| ~ ^ units of this asset 
for its price x. A unit of this asset has an expected value D some time later. Hence at this time, 
the money m invested in asset A has value 

M = — m 

X 

If the agent is lucky, then at the expiration time — > 1, and his money m has become valuable by 
a factor 

_ M _ D 
^ m X 

It is reasonable to assume that the agent wants to spend his money in an asset of which he expects 
that > 1. Thus, depending on his expectation about the future value of the agents buy or 
sells this asset. Therefore, if the agent expects that > 1, he will buy, otherwise he will sell. This 
causes an increase (decrease) of demand in this asset. Due to the increase (decrease) of demand, 
the price will rise (fall). Therefore, price evolution is thought to depend on the expected growth 
rate. 

Theoretical results due to our model are compared with empirical data from daily returns in 
the period from 01/01/1990 to 112/31/2004 of 10 large indices listed in table [l] Resuffs concern 
i.) the non-gaussian distribution of log-returns, ii) the non-linear Hurst spectrum of return trials, 
and iii.) the Renyi entropy. 
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France 


CAC 40 


2 


Germany 


DAX 30 


3 


Hong Kong 


Hang Seng 


4 


Japan 


NIKKEI 225 


5 


Switzerland 


Swiss Market Index 


6 


Switzerland 


Swiss Performance Index 


7 


United Kingdom 


FTSE 100 


8 


United States 


Dow Jones IndAvg 


9 


United States 


Nasdaq 100 


10 


United States 


S&P500 



Table 1: For each index, we considered daily data from 01/01/1990 to 12/31/2004 provided by 
Thompson Datastream. 
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1 An expectation driven market 



The price process X of an index from time to time t is described a the concatenation of n 
independent trading periods = [r, r + 1), where '[r/ and ',t + 1)' can be regarded as the 
opening time or the closing time of this period. The process from then is 

[0, t) = [0, 1) 2) t) 

Let X be the closing price of period Tt and x' the opening price of the next period. Are x and x' 
independent from each other? In reality there is night between both periods in which a lot can 
happen. On the other hand, reality tells us that the opening price does not differ too much from 
the former closing price. Now, we assume that the process can be described as a concatenation of 
independent processes defined on periods T-^. This means that we assume that prices x = X^) and 
x' = X\^tj^i are independent. 

We consider a simple 1-period model of a financial market with only one asset, whose 

price at time t \s x > 0. Prices on this market are supposed to follow a multiplicative random 
process given by 

x' = Tx (1) 

where x' is the price at time t' . The growth rate F > is assumed to be due to the expectations 
investors have at this time about future prices: Investors are supposed to build their believes about 
the growth rate based on some economic entity available today. A particular simple setting is that 
the expected growth rate is an increasing function of an expected dividend yield or earnings rate 
y, i.e. 

r = r-(.) :=$(f^ 

while $(0) = 0. D represents is a non-negative random variable which is distributed according to 
some (stationary) distribution in some finite interval [0, d\. 

In the following we assume that for all $ is a power law, i.e. 

K?) - m 

where a > is a constant. The second assumption is about the distribution J^: 

D^U{Q,d). (3) 

This assumption is due to the lack of knowledge about the entity D. As well known, if all values 

within the interval are taken with equal probability, the uniform distribution is the unique distribu- 
tion that minimizes information. This can be taken as the statement that all possible information 
are arbitraged away from a financial market. 
With 

D = Sd, 6€U{0,1) 

and y{x) := ^, we finally obtain 

x' = X ^6 y{x) 

Given x, the probability that the gross return i? = ^ is larger than some r is given by 

FR{r\x) = P[R > r\x] = P[(^y(a;))" > r \ x] 

y{x) 
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for < r < (^)"- Un-conditioning by integrating over all (initial values) x implies that the 
unconditioned cumulative tail distribution obeys a power law 



Mr) = 




Therefore log-returns Z — In R are double exponentially distributed according to 

fziz)^:^e~i\^\. (5) 
z a 

For further improvement of the model, it might be interesting to list the assumptions made. 

1. The financial market contains only 1 asset; 

2. Due to equation 21 the result holds for a 1-period model; 

3. Due to equation 13 payoffs are uniformly distributed within a fixed finite interval; 

4. Due to equation 121 the growth rate is r(a;) is a power law with a constant scaling exponent 
a. 

Hence the model must be regarded as a 0-ORDER APPROXIMATIONS. Particularly, a financial 
market contains more than one asset, see (1), the price evolution is a process rather than can be 
described by a series of independent periods, see (2), (3) is reasonable only if one assumes that for 
all times all values in [0, d] are equally probable. 

Assumption (4) has an interpretation in terms of market liquidity. Let us assume that the 
volume traded is a function of the price 

V = V{x) ~ - 

X 

i.e. the more costly the asset is the smaller is the amount of this asset that can be bought by one 
unit of capital. Then the price change Z — Inx' /x generated by a trade of volume V{x) at price x 
yields 

x' 

Z{x) = In — = In 
~ a hiV{x), 

up to some additive constant. Therefore the price impact of a trade of size V at price x equals 

dZ{x) = ad\nV{x). (6) 
Hence if a w 0, a large trading volume is necessary to move the price. 

Therefore a can be regarded as a kind of elasticity of the market at the price level x, that 
might economically understood as a measure of current market depth or its liquidity 0. In the 
light of this interpretation of the parameter a, the fourth assumption means that liquidity should 
be independent of the price level and moreover constant for any time period. 

Being aware of these assumptions, we compare our result, equationEl with real data. Deviations 
of empirical return distributions from our theoretical result shall indicate which assumptions we 
have to modify to get a better model. 
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2 Numerical estimates of "old friends" 

Table n displays the set of indices considered"'^. Dealing with indices, makes it necessary to discuss 
the meaning of the entity D. Originally D was meant as representing the process of expected 
dividend pay off. Indices do not pay off dividends. Staying close to this interpretation, D has 
to be regarded here as a sum of expected dividend flows of the constituting assets. Furthermore, 
expected dividend are not the only source of signals that affect expectations about the growth 
rate of an index. Thus signals are usually summarized as 'news' in the economic literature. Hence 
D could also be regarded as representing the news process in the market. The broader picture 
therefore is that D represents the process of any signal that affects expectations about growth 
rates of the respective market. The first question to be considered is to which degree our theo- 
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Figure 1: Returns distributions of the S&P 500, left: a semi-logarithmic plot, right: a double logarithmic 
plot 

retical result, equation agrees with empirical data. Let us first have a look at the distributions 
at the S&P 500, see Figure Q] The left picture is a semi-logarithmic plot and the right picture 
is a double logarithmic plot of daily returns from 01/01/1990 to 12/31/2004. Cross hairs mark 
negative returns and dots display positive returns. Recall that an exponential distribution in a 
semi-logarithmic plot becomes a straight line, while a straight line in a double logarithmic plot 
corresponds to a power-law distribution. Numerical estimates of the parameter can be found in 
Appendix ^ 

Quantile-Quantile plots provide a good descriptive method to judge about whether two sets of 
data come from the same distribution. A line indicates that data in both sets are very likely to 
come from the same distribution. As an example, let us consider the S&P 500, see Figure 13 

Figure [3 shows the QQ-plot of the empirical time series of log returns Z with respect to an 
exponential distribution. The left column considers positive log-returns := ZyQ, while the 
right column is for negative log-returns Z_ Z<q. The first row displays the QQ - plots of Z± 
which respect to the exponential, while the second row displays a QQ-plots for In \ Z± \ with respect 
to the logarithm of the exponential. The scheme see in Figure El thus is: 
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In this version, only the results for S&P 500 are shown. They are in fact representative for the others considered. 
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QQplot positive returns 
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Figure 2: QQ-plots of the S&P500 with respect to the exponential distribution 



The QQ plots of all indices considered show the same pattern: While the middle part of each 
QQ-plot is linear, systematic deviations from the straight line in the QQ-plots occur for either 
small returns or large returns. Particularly, deviations from the diagonal for the QQ-plot of Z± 
are seen for high quantiles, while deviations form the line in the QQ-plot of In \ Z± \ exists for small 
quantiles. This means that empirical returns deviate from being exponentially distributed for very 
large returns, see the first row, and for very small returns, see the second row. This agrees with the 
observations in the respective pdf 's, where we see that typically the empirical pdf has less mass in 
than the Laplacian, while it's tails are usually fatter than those from the Laplacian. 

Fgures show the degree to which our elementary mouse model fits to empirical data: Deviations 
from a clear linear relation exist either for small returns and for large returns, while in the middle 
range, this method indicates that the proposed double exponential distribution provides a fairly 
good description of the empirical data. 
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3 On multiscaling in time series of log-returns 



Stylized facts are important statistical properties since they are seen in the empirical returns of 
'all' financial markets. Besides the non-Gaussian character of empirical asset returns distributions, 
a second important fact concerns multiscaling in the time series of empirical asset returns. The 
existence of multiple scales in the system implies that returns distributions are not invariant under 
the choice of different time-scales, i.e. one observes that the distribution of returns with lags of 
the order of minutes, days, week, and so forth deviate from each other, see IH]. 

In the following we estimate the scaling exponent ('Hurst exponent') 



m\zit,T)\' 



from time series of returns Z{t, T) = In ^^^^^ with time scale T. In our 1-period model, T — 1. It 
can be shown that r(0) = — 1, while r is increasing and concave. For a simple (fractional) diffusion 
process, r is a linear function of q. The deviation of r from being linear, is an important issue in 
determining the multi scaling nature of the underlying process. Figure |31 show T{q) as a function 
of q for various indices compared with our model. For orientation, the dotted straight line has a 
slope to 1/2, which corresponds to Gaussian diffusion, while the solid line is the estimated graph 
of scaling exponents in our model. Multi-fractal processes have been proposed as a new formalism 



Hurst spectrum MODEL 



Hurst spectrum DATA 




02468 10 02468 10 

Figure 3: non-linear HurSt spectnmr in the time series of daily returns of ttft S&P500 (right), to 
be compared with our model with the estimated a, see tables HJ O (left) 

for modeling the time series of returns in finance. The major attraction of these processes is their 
capability of generating various degrees of long-memory in different powers of returns - a feature 
that has been found to characterize virtually all financial prices, see [lUllTIIT]. The prominent issue 
of these modeling approaches is the use of cascading processing. 

The multifractal formalism is interesting even from another point of view. It may open the 
door wider to bring two disciplines closer to each other: Statistical Physics and Financial Markets, 
the connection being the Reny entropy'. For an introduction into the field of thermodynamics and 
non- linear systems and related concepts see [2]. 
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4 The Renyi Entropy of our stylized financial market 



Time series of log-returns show multiscaling. The Renyi entropy has proven to be a reasonable 
entropy measure for multifractal systems in which long-range correlations exists. We therefore 
estimate the Renyi entropy defined by 



1-/3 ^ 



(7) 



for a system with r micro states and a zooming parameter /3 G M. Figure^ shows the Renyi en- 
tropies for different indices. We will compare the theoretical Renyi entropy with the Renyi entropy 

3.5 



2.5 



1.5 



DAX 30 
DJIA 

HangSeng 

NASDAG 

NIKKEI 

SP500 

SMI 

CAC 40 
FTSE 



0.5 

01 23456789 10 

Figure 4: The Renyi entropies of single indices as a function of f3. 

found in empirical data. 

We consider a trail of log-returns Z and partition the range into r > cells each of length 1/r. 
Then the k-th cells is 

"fc- 1 k] 



The probability pk that Z € Ck equals 



k ~ l..r. 



Pk 



fz{x) dx 
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According to equation [SJ Z is Laplacian distributed. Therefore we obtain 

h. 

1 



Vk = 



dx 



Summing probabilities over the sample space {1, .., r} gives 



dx — — 



'dx 



1 - e" 



For later purposes we define 



1 — e a 



(8) 



Therefore we normalize probabilities by 



TTfc = Ci e 



The Renyi entropy yields 

Mp) 



1-/3 ^ 1-13 ^ ^ 

k 

InCf + ln^e-° 



1 



From 2^^^ e 



k 



1-/3 

we obtain for the Renyi entropy 



(9) 



To estimate the parameter a from measuring the Renyi entropy, it is therefore sufficient to consider 
two cases, /? = 0, 1, where -Ro(p) ■= hm^^o Rpip) ^^id Ri{p) := lim^^i R0{p). From its definition 
it immediately follows the well known relation 



Ro{p) = Inr 



(10) 



Further, in the limit /3 — > 1, the Renyi entropy becomes the Shannon-Boltzmann entropy: Ri{p) 
— TT; InTTfc and therefore 

^i(p) ^ -^Cie~^ ln(Cie- = 

k 



InCi 

a r 



ar\2 



so that we obtain 



.lnCi + -^y k = -lnCi + ^^'' 

a T — ^ rv r 

k 



R,{p) = -InCi + ^Ci. 

2 a 



(11) 
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Since Ci = Ci(a,r), equations 1101 and 1111 allow to estimate the distribution parameter a from 
empirical data - under the hypothesis that the distribution of Z is exponential. 



We calculate the Renyi entropy for return trails of each single index 1, .., 10, see Table for a 
fixed number of cells r = 30, see Figure|Sl Former estimates of the asymmetry of the distribution of 
negative returns and positive returns, measured by the parameters a± showed that the distribution 
is (almost) symmetric. Hence we considered the trail of absolute returns, i.e. | lnZ|. We normalized 
returns to the unit interval by considering \Z*\ ~ max izY — ■ l^ft upper picture shows the 

distribution of \Z*\ and a fit with respect to the exponential distribution given the (adjusted) a 
estimated from the exponential distribution as in Appendix^ The upper right picture shows the 
resulting empirical Renyi entropy (solid line) and the graph of Rp (p) from equation|51 (dashed line) . 
The prediction of our model agrees fairly well with the data. The logarithmic plots in the lower 
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Figure 5: 'Sii:P 500: Distribution of absolute returns and Reny'fe&tropy 



row show the deviations for small j3 G [0, 3] and large (3 e [25, 50], see Figure El Recalling that P 
is a zooming factor, which places emphasis on different probability regimes, these deviations agree 
with the findings stated earlier: Our model gives too much mass to high probabilities, related to 
small returns, while it does not give enough mass in the tails, i.e. for large returns. 
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5 Conclusion 



All models are wrong but some are useful. 

- G. E. P. Box, 1979 

Our aim is to understand price dynamics on a financial market. The existence of stylized facts 
suggests that price trails of different financial markets might be regarded as different realizations of 
a more general stochastic system, called 'The financial market'. If so then the question is about the 
nature of this system. Since prices are macro-observables of a financial market, the model about 
price dynamics is defined on the macro level. Due to the set of assumptions used in its derivation, 
this model is an approximation in itself. Results therefore can also only as zero-approximations as 
indicated by the list of assumptions made. We estimated three major properties: the distribution 
of (logarithmic) asset returns, the Hurst exponent in their time series, and finally the Reny entropy. 
Although the model is a zero-order approximation, theoretical results are already in fairly good 
agreement with real data. First-order corrections concerning the set of assumptions made should 
improve these theoretical findings. Taking all this together, this model might serve as a good 
starting point for further improvements. These steps should follow from observing where our 
theoretical results deviate from empirical data and by successively and modestly modifying the 
assumptions made. 
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A Quantitative estimation of parameters 



We compare two distributions with respect to their goodness to fit empirical data. One is the 
exponential with parameter a, while the other one is the WeibuU with parameters (a, fj,) 

a 



/rw = -(-) 

z V a/ 



The exponential distribution is obtains from the WeibuU in the case ^ = 1, while in the limit /i ^ 0, 
fz{z) approximates the Pareto distribution arbitrarily well, see jHj. Fits of the distributions for 
positive and negative returns due to these distributions therefore give respective parameters 

/3 = (a^,m^ ; a~,m~) . 
To make differences more obvious, we standardize log returns in the usual way 

^, _ Z± - {Z±) 

where Z'^ are the positive and negative returns, respectively, and a denotes the respective standard 
deviation. Therefore the estimates parameter /i is of special interest. In the following we summa- 
rize log likely fits of positive and negative returns, respectively, to the exponential distribution and 
the WeibuU distribution, giving parameters and (a^, respectively. 

Distributions are close to being symmetric a+ ~ q;_, see tables |21 and |5| Furthermore ^ 
« 1, i.e. their distributions are close to an exponential distribution. This finding is also 
supported by considering the respective entropic distances between the empirical distribution and 
the exponential and WeibuU distribution respectively as formalized by the Akaike's Information 
Criterion 

AIC ^-2 log {c{9\x) \ +2K, 



where K is the number of parameters, here K=l for the exponential and K = 2 ioi the WeibuU 
distribution. The term 2 K can be considered as a penalty for introducing additional parameters. 
C{9\x) is the maximum log-likelihood of the parameter 9 given the data x. AIC is in particular 
useful for nested model such as the exponential and the WeibuU in this case. However, we considered 
the ratio of the AIC's of both models 

C{exponential) + 1 
^^''^ " C{Weibull) + 2 ' 

As seen in tables O and |3 the ratio is close to 2/3, which indicates that the approximation that 
the distribution of log returns is an exponential is quite well and that we can trust the estimated 
parameters ak- 
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INDEX 


#Z<0 


a 


a 




ratio 


DAX 30 


1806 


0.6985 


0.6613 


0.8885 


0.6671 


SWISS SMI EXP 


1066 


0.7021 


0.6762 


0.9209 


0.6673 


FRANCE CAC 40 


1849 


0.7054 


0.6869 


0.9399 


0.6671 


FTSE 100 


1840 


0.7074 


0.6997 


0.9742 


0.6671 


SWISS SPI EXTRA 


955 


0.7342 


0.6891 


0.8820 


0.6673 


DOW JONES 


1793 


0.6891 


0.6713 


0.9437 


0.6671 


Hang Seng 


1808 


0.6393 


0.6000 


0.8815 


0.6671 


NASDAQ 100 


1745 


0.7077 


0.6887 


0.9405 


0.6671 


NIKKEI 500 


1864 


0.7500 


0.7698 


1.0688 


0.6671 


S&P 500 


1799 


0.6869 


0.6663 


0.9353 


0.6671 



Table 2: Parameter estimation for negative returns of the indices considered 



INDEX 


#Z>0 


a+ 


a+ 




ratio 


DAX 30 


1964 


0.7057 


0.7280 


1.0840 


0.6671 


SWISS SMI EXP 


1190 


0.7003 


0.7186 


1.0683 


0.6673 


FRANCE CAC 40 


1925 


0.7454 


0.7786 


1.1285 


0.6671 


FTSE 100 


1948 


0.7339 


0.7638 


1.1129 


0.6671 


SWISS SPI EXTRA 


1283 


0.6634 


0.6984 


1.1615 


0.6673 


DOW JONES 


1979 


0.7223 


0.7478 


1.0948 


0.6671 


Hang Seng 


1904 


0.7038 


0.7105 


1.0229 


0.6671 


NASDAQ 100 


2033 


0.6933 


0.7070 


1.0485 


0.6671 


NIKKEI 500 


1825 


0.6592 


0.6173 


0.8699 


0.6670 


S&P 500 


1983 


0.7207 


0.7373 


1.0589 


0.6671 



Table 3: Parameter estimation for positive returns of the indices considered 
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